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X-RAY LOCALIZATION OF THE GLOBULAR CLUSTER Gl WITH XMM-NEWTON 
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ABSTRACT 

We present an accurate X-ray position of the massive globular cluster Gl by using XMM-Newton and the Hubble 
Space Telescope (HST). The X-ray emission of Gl has been detected recently with XMM-Newton. There are two 
possibilities for the origin of the X-ray emission. It can be either due to accretion of the central intermediate-mass 
black hole, or by ordinary low-mass X-ray binaries. The precise location of the X-ray emission might distinguish 
between these two scenarios. By refining the astrometry of the XMM-Newton and HST data, we reduced the 
XMM-Newton error circle to 1.5". Despite the smaller error circle, the precision is not sufficient to distinguish an 
intermediate-mass black hole and luminous low-mass X-ray binaries. This result, however, suggests that future 
Chandra observations may reveal the origin of the X-ray emission. 

Subject headings: binaries: close — globular clusters: individual (Mayall II = Gl) — X-rays: binaries 



1. INTRODUCTION 

Globular clusters are very efficient places to produce X-ray 
binaries via dynamical interactions. In particular, it has been 
known for many years that the formation rate per unit mass of 
luminous (Lx > 10 36 ergs s" 1 ) X-ray sources is much higher in 
globular clusters than that of the rest of our Galaxy. More re- 
cently, similar results are found in other nearby spiral galaxies 
like M31 (DiStefano et al. 2002) and M104 (DiStefano et al. 
2003). Among all extragalactic globular clusters, Gl in M31 is 
an intriguing one. With a luminosity of ~ 10 6 L Q (Rich et al. 
1996), it is the most luminous star cluster in the Local Group, 
and also one of the most massive, at (7- 17) x 10 6 M Q (Mey- 
lan et al. 2001). The rates at which X-ray binaries are cre- 
ated in the cluster core are therefore expected to be high com- 
pared with globular clusters in the Milky Way. Furthermore, 
it has been claimed, based on kinematic studies, that Gl hosts 
a~2x 10 4 M Q intermediate-mass black hole (Gebhardt et al. 
2002,2005). However, this result is controversial and has been 
challenged by Baumgardt et al. (2003). X-ray observations 
of Gl therefore allow us to investigate some of the interesting 
properties of the cluster. 

Recently, XMM-Newton has conducted three short (< 10 
ksec) observations of Gl and has discovered an X-ray source 
coincident with Gl (Trudolyubov & Priedhorsky 2004; Pooley 
& Rappaport 2006). To explain the origin of the X-ray emission 
of Gl, Pooley & Rappaport (2006) proposed that it could be 
due to accretion of ionized cluster gas by a central intermediate- 
mass black hole or it could be produced by a conventional X-ray 
binary, and it is also possible to distinguish these two scenarios 
by obtaining a precise localization of the X-ray emission. If 
the X-ray emission is due to a central 20,000 M© black hole, 
we expect it comes from within 50 mas of the center. How- 
ever, if low-mass X-ray binaries are responsible to the X-ray 
emission, then we expect the X-ray emission offsets from the 
core. This requires high-resolution X-ray observations. How- 
ever, there is no Chandra observation of Gl and only XMM- 
Newton observations are available. Although Pooley & Rap- 
paport (2006) investigated the XMM-Newton spectra of Gl in 
detail, they did not perform an astrometric study. The absolute 
astrometry of XMM-Newton is about 2" (Kirsch 2006) while 
the statistical uncertainty is intensity dependent. This leads to a 



positional error of about 2" -6" depending on the source bright- 
ness. While the spatial resolution of XMM-Newton is much 
poorer than Chandra it is possible to localize X-ray positions 
to I" — 2" with XMM-Newton if one calibrates the astrometry 
carefully. In this paper, we refined the X-ray position of Gl by 
performing precise relative astrometry using XMM-Newton and 
the Hubble Space Telescope (HST). 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. XMM-Newton 

Gl was first observed with the XMM-Newton in 2001 Jan- 
uary for a total exposure time of ~ 8 ksec. There were two more 
XMM-Newton observations in 2002 December and 2003 Febru- 
ary. Both observations were off-axis resulting heavy vignetting 
and one of the observations was affected by high background 
[see Pooley & Rappaport (2006) for a summary]. In this Let- 
ter, we only consider the first observation taken in 2001. All 
three cameras (one pn camera and two MOS cameras) of the 
European Photon Imaging Camera (EPIC) were turned on for 
collecting data. All the X-ray data were processed with the 
XMM-Newton Science Analysis System (SAS) version 7.0. 

We downloaded the raw data from the XMM-Newton archive 
and reprocessed with SAS together with the latest calibration 
products. The reprocessed event lists were first examined for 
background variation using the high energy (10-15 keV) back- 
ground lightcurves and we did not find any significant back- 
ground flaring event. We extracted X-ray images with pho- 
ton energies in the range of 0.3-10 keV, and only considered 
events with FLAG = and single and double events for the pn 
camera (PATTERN < 4), and single to quadruple events for 
the MOS cameras (PATTERN < 12). Source detection was 
then performed using a maximum likelihood approach as im- 
plemented by the SAS tools edetect_chain. We ran the source 
detection simultaneously on the data from all three cameras. Gl 
was clearly detected and was seen in all three cameras with a 
combined detection likelihood of 63. We compared the X-ray 
source list with the 2MASS and USNO catalogs and images, 
and looked for coincidence of bright and isolated stellar ob- 
jects. We found one star (2MASS 00325251+3931424) that is 
< 3" from the X-ray position and it is likely to be a foreground 
star. 
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To verify the nature of this X-ray emitting stellar object, we 
computed the hardness ratios. These ratios were based on the 
source counts in the three energy bands: S (0.3-1 keV), M (1-2 
keV), and H (2-10 keV). The two hardness ratios are defined as 
URl=(M-S)/(M+S)andUR2=(H-S)/(H+S). Figure 1 shows 
the color-color diagram of the X-ray emitting foreground star 
and Gl. We have overlaid the color-color diagram with lines 
showing the tracks followed by representative spectra with dif- 
fering values of Nh- The X-ray colors of the X-ray emitting star 
indicate that it has a very soft X-ray spectrum, consistent with 
a very soft X-ray source (Di Stefano & Kong 2004). The X-ray 
radiation is therefore likely due to the coronal emission from 
a foreground star. The star has a R magnitude of 14.0 (Monet 
et al. 2003). We calculated the X-ray to optical flux ratio as 
log(/x//*) = log/ x + 5.67 + 0.4/? (Hornschemeier et al. 2001). 
With a count rate of 0.018 c/s in the pn detector and assuming 
a Raymond-Smith model with kT RS = 0.3 keV and Nh = 10 21 
cm" 2 , the 0.3-10 keV flux is 2.8 x 10" 14 ergscm" 2 s _1 and the 
corresponding fx /fit is 0.005, typical for a foreground star 
(Hornschemeier et al. 2001). 

Based on the optical counterpart, the boresight correction 
that needs to be applied to the X-ray source positions is 1 .47" ± 
0.86" in R.A. and 2.51" ±0.86" in decl.; the uncertainties are 
quadratic sum of the positional errors of the X-ray and 2MASS 
source. The correction is consistent with the absolute point- 
ing accuracy of XMM-Newton (see also an example in Pietsch, 
Freyberg & Haberl 2005). 
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FIG. 1 . — Color-color diagram of Gl and a nearby X-ray emitting object 
(lower left). Also plotted are the estimated hardness ratios estimated from dif- 
ferent spectral models. Top to bottom: Power-law model with a of 1.0, 1.5, 
2.0, and 2.5 and Raymond-Smith model with kT R $ of 0.3 keV. For each model, 
Nh varies from the left from 5 x 10 20 , 10 21 , and 5 x 10 21 cnr 2 . 

2.2. Hubble Space Telescope 

Gl was observed with the HST Advanced Camera for Sur- 
veys (ACS) in High Resolution Channel (HRC) mode on 2003 
October 24. The total integration time is 41 minutes in the 
F555W filter centering on Gl. We used the HST pipeline data 
that were shifted and co-added using the MultiDrizzle package 
in PyRAF, with masking of cosmic rays, saturated pixels, and 



bad pixels. We calibrated the astrometry of the HST data by 
using the 2MASS catalog. The field-of-view of the ACS/HRC 
is small (29" x 26") and only two stars in the field are in the 
2MASS catalog. By computing the average offset between the 
HST and 2MASS stars, we shifted the HST image by 0.77634" 
in R.A. and -0.5814" in decl., with a residual of 0.13". The 
ACS/HRC image of Gl is shown in Figure 2. 

Gl was also observed with the HST Wide Field Planetary 
Camera 2 (WFPC2) on 1995 October 2 with a total integra- 
tion time of 37 minutes in the F555W filter. In addition 
to the F555W data, images were also taken in the F814W 
and F1042M filters. We downloaded the F555W image from 
the WFPC2 Associations Science Products Pipeline for which 
cosmic-ray free, science-quality images are dithered and co- 
added. Since the field-of-view of WFPC2 is much larger than 
that of ACS/HRC, we can correct the astrometry with 6 stars in 
the 2MASS catalog. We applied the astrometry correction us- 
ing IRAF task ccmap yielding a residual of 0.14". The WFPC2 
F555W image is shown in Figure 2. 

3. X-RAY LOCALIZATION OF Gl 

After registering the absolute reference frames of the HST 
and XMM-Newton images to the 2MASS catalog, we located 
the center of Gl in the HST images and the X-ray position of 
Gl. We determined the centroid of Gl in the ACS/HRC im- 
age (R.A.=00h32m46.537s, decl.=+39d34m40.65s with la er- 
ror of 0.004") by computing the intensity weighted mean within 
the core radius (0.21"; Ma et al. 2007) using IRAF task cen- 
ter. We also checked the result by using the half-mass radius 
(1.73") and there is no difference except for a larger error bar 
(0.01"). For the X-ray position, we applied the astrometric 
correction on the value determined by edetect_chain yielding 
R.A =00h32m46.6s, decl.=+39d34m40s. We then determined 
the la radius error circle (1.5") of the XMM-Newton position 
of Gl by computing the quadratic sum of the positional uncer- 
tainty for the X-ray source ( 1 .23"), the uncertainty in the optical 
astrometry (0.13"), and the uncertainty in the X-ray boresight 
correction (0.86"). The same procedure was also applied to the 
WFPC2 image. Figure 2 shows the ACS/HRC and WFPC2 im- 
ages of Gl and the la radius X-ray error circles. The center of 
Gl derived from the optical image is also marked. 

4. DISCUSSION 

Using XMM-Newton and HST, we determined the centroid 
of Gl in the optical images as well as the X-ray position of Gl. 
From Figure 2, although the X-ray position offsets from the 
cluster core, the cluster center is within the 1 a error circle of the 
X-ray position. Therefore, the current XMM-Newton data can- 
not provide constraint on whether the X-rays of Gl come from 
Bondi accretion of ionized cluster gas by a central intermediate- 
mass black hole for which the X-rays should come from the 
central 50 mas of the cluster (Pooley & Rappaport 2006). Al- 
ternatively, the X-ray emission could be produced by luminous 
low-mass X-ray binaries and we expect such X-ray emission 
may be outside the cluster center. Previous X-ray observations 
of globular clusters suggest that luminous (L x <; 10 36 ergs s" 1 ) 
X-ray sources tend to locate within the core radius (Grindlay et 
al. 1984). Recent Chandra observations also show that nearly 
half of the quiescent low-mass X-ray binaries are found within 
the core radius (e.g., Grindlay et al. 2002; Pooley et al. 2002; 
Heinke et al. 2003,2006). Therefore, it is likely that a lumi- 
nous low-mass X-ray binary would locate within the core ra- 
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FIG. 2.— HST ACS/HRC F555W (left) and WFPC2 F555W (right) images of Gl. The cluster core is marked by across. The circles of both images are the la 
radius (1.5") XMM-Newton error circles. The two bright stars in the field were used for calibrating the astrometry with the 2MASS catalog. 



dius (0.21"; Ma et al. 2007) of Gl. The X-ray emission of Gl 
could also come from multiple low-mass X-ray binaries and 
we may be able to resolve Gl as an extended source with high 
spatial resolution instrument. Pooley & Rappaport (2006) esti- 
mated that about 75 low-mass X-ray binaries might be in Gl. 
It is worth noting that only one globular cluster, Ml 5, is known 
to host two luminous X-ray sources (White & Angelini 2001; 
Hannikainen et al. 2005). In conclusion, based on the cur- 
rent XMM-Newton data, we cannot distinguish the two possible 
mechanisms of generating the X-ray emission of Gl. 

While the X-ray position of Gl is the most crucial factor to 
determine its nature, Pooley & Rappaport (2006) also suggested 
that the X-ray spectrum may provide some hint. However, it has 
been proven that X-ray spectra of intermediate-mass black hole 
candidates consist of a class of different spectral shapes and in 
many cases, the spectra can be fit with several models and the 
estimated mass of the accreting black hole is model dependent 
(e.g. Stobbart et al. 2006; Goncalves & Soria 2006). Some 
intermediate-mass black hole candidates also show X-ray spec- 
tral change at different luminosity states (Kong & Di Stefano 
2005). These make interpretation based on X-ray spectra more 
difficult. Nevertheless, if we assume a simple accretion disk 
model for Gl, following Pooley & Rappaport (2006), we would 
expect Gl has a 10 eV supersoft component (see Di Stefano & 
Kong 2003); any emission above 0.5 keV must come from ad- 
ditional components. From the color-color diagram (Figure 1), 
Gl has significant emission above 1 keV and indeed it is very 
similar to a typical X-ray binary in M3 1 with a simple power- 
law spectral model (Kong et al. 2002). Therefore, if Gl has a 
10 eV supersoft spectrum, it must also have an additional hard 
component. Indeed, it would be a challenge for XMM-Newton 
and Chandra to detect such supersoft emission because an ab- 
sorbed (Galactic value to the direction of M31; N H = 7 x 10 20 
cm -2 ) 10 eV spectrum turns over at about 0.2 keV which is the 
sensitivity limit of these instruments. For instance, if the X-ray 
emission is dominated by a 10 eV spectrum, simulation shows 
that it requires 700 ksec XMM-Newton or 1 Msec Chandra ob- 
serving time in order to detect the source. If the black hole of 
Gl is only 100 M Q , the thermal emission would have a tem- 
perature of about 80 eV and we should be able to detect it with 



a 6 ksec XMM-Newton or 10 ksec Chandra observaiton. As a 
supersoft X-ray source (kT < 100 eV), we do not expect to see 
X-rays above 1 keV which is not consistent with our current 
result. Alternatively, if the X-ray emission is from a luminous 
low-mass X-ray binary, we would also expect soft multi-color 
disk blackbody X-ray emission (kT in w 0.3 - 3 keV) in addition 
to a power-law like component associated with Comptonization 
of cooler photons (e.g. Sidoli et al. 2001). However, with only 
~ 70 counts from all three XMM-Newton detectors, we do not 
have a good constraint on the spectral model. In conclusion, 
X-ray spectra provided by XMM-Newton and Chandra alone 
cannot provide convincing evidence for the nature of the X-ray 
emission from Gl. 

Although the X-ray position provided by XMM-Newton can- 
not provide any reasonable constraint to the nature of the X-ray 
emission of Gl, it suggests that future Chandra observations 
may resolve the problem. As discussed in Pooley & Rappa- 
port (2006), we can improve the relative astrometry of Chan- 
dra and HST to 0.1" -0.2" if we can match a few Chandra 
sources to their optical counterparts (e.g. foreground stars or 
background active galactic nuclei). With such observations, we 
can accurately localize the X-ray emission of Gl. However, the 
fact that a luminous low-mass X-ray binary is likely to locate 
< 0.21" from the cluster core suggests that we may not disen- 
tangle from the emission of a possible intermediate-mass black 
hole. Alternatively, a Chandra observation may be able to dis- 
tinguish between multiple X-ray sources as an extended object 
and point-like emission. 
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